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Abstract 

The time-evolution of quantum meson fields in 0(4)-linear sigma model is treated 
approximately. It is shown that the amplification of the amplitudes of pion modes 
with low momenta occurs by means of both the parametric resonance and the forced 
oscillation. 
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It is interesting to study the dynamics of the chiral phase transition in the context of 
relativistic heavy-ion colhsions, such as the problem of the formation of a disoriented chiral 
condensate. One of the important theoretical aspects is to investigate the time-evolution 
of the order parameter of the chiral phase transition. As for the fluctuation modes around 
the order parameter, the amplification of the amplitudes must occur accompanying with 
the relaxation of the chiral order parameter. Recently, the present author has investigated 
the time-evolution of a collective meson field in the context of the dynamical chiral phase 
transition. It has been seen that the amplitudes of quantum fiuctuations with low momenta 
around the mean field configuration have been amplified with the direction of pion modes 
in the 0(4)-linear sigma model. As is indicated by many authors, this phenomena may 
be understood in terms of a parametric amplification.^)"^) However, there may be another 
mechanism to amplify the fiuctuation modes. Actually, it is seen that the amplitude of the 
pion mode with higher momentum without k = can be amplified in the previous paper. ^) 
In this paper, it is pointed out that there is another mechanism to amplify the quantum 
fluctuation modes through the chiral phase transition, namely, that a forced oscillation 
works as well as a parametric resonance. 

Before treating the 0(4)-linear sigma model, it is instructive to recapitulate the ingre- 
dients of the parametric resonance and the forced oscillation. Let us consider the following 
equation of motion for x{t) : 

X + u;l{l — hcos^t)x — . (1) 

The parametric resonance occurs around 7 = 2uJo/m (m = 1,2, •••).^) For m — 1, the 
resonance region of the frequency 7 is given as 

-huj/2 < e < huj/2 , for 7 = 2a;o + e , e < , (2) 

where we treated h as a small parameter, h <^ 1. Similarly, for m = 2, we obtain 

-(5/24) • /i^c^o < e < (1/24) • /i^c^o , for 7 = u;o + e , e < a;o ■ (3) 

In general, for given integer m, the unstable solution of x{t) is obtained as the form x{t) oc 
e** X (oscillation and constant parts), where s is of the order of h'^. Namely, the amplification 
becomes slower in time as the integer m becomes larger. The resonance region also becomes 
narrow with the order of k^, that is, —0{h^) < e < 0{h"^). If h is not so small, we have to 
deal with the equation (1) directly. The equation (1) is known as Mathieu's equation and 
the properties of this equation are investigated in detail. In the following investigation, h 
is actually small fortunately. Thus, it is not necessary to take care of the full treatment of 



2 



this Mathieu's equation. Next, let us consider the equation of motion which describes the 
forced oscillation : 

X + ujIx = F cos . (4) 

The solution with initial conditions x{0) = i;(0) = is obtained as 

,, ^ cos i7i — cos a;oi 2F . ouq — f2 . u!o + f2 

xit) = h ■ ^5 -— = -. — r-^ — r- sm — - — t ■ sm — - — t . (5) 

As is well known, when Q has a value near uq, the behavior of the above solution shows the 
beat. Especially, i? « a;o, the amplitude of the solution becomes large linearly in time : 

x{t) {F/2ujq) ■ tsvaujot (6) 

Now let us return to the main subject, namely, let us treat the 0(4)-linear sigma model in 
the squeezed state approach. We start from the following Hamiltonian density of 0(4)-linear 
sigma model: 

V. = \7ta{xf + \V(t>a{x) " V(t)a{x) + A ((t>a{xf - / - Hct>,{x) . (7) 

The squeezed state is adopted as a trial state in the time-dependent variational principle : 

\<P{t)) = Arexp(i(7f •(/.-^•7r))exp(0- [-(G-^ -G(°)-^)/4 + ir] -0)10) , (8) 

where A'^ represents a normalization factor and W ■ cf) = J2a=o I d^Wa{x, t)(t)a{x) and • G^^ ■ 
= T.l=ond?xdhj4>a{x)G-\x,y,t)4>a{y) and so on. Here, Of = (O|(/.„(x)0„(2/)|O). The 
expectation values of the field operators are obtained, for example, as 

{^{t)\Mx)My)Wit))-^aix,tpaiy,t) + Ga{x,y,t) . (9) 

^^(a;, y) and Ga{x, y, t) represent the mean field and the fluctuation around it, respectively. 
The time-dependence of these functions, together with Tf^ and E^, are determined in the 
time-dependent variational principle : 



6 j^ymm- - H\<p{t)) = , (10) 



where H = ^SxTi. The two-point function Ga{x,y,t) can be further expanded in terms of 
the mode functions under the assumption of the translational invariance : 



Ga(x,y,t)^ jj^^^^-y^litf , (11) 
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where = / (i'^k/(27r)^. The equations of motion for the condensate iPa{t) and the fluctua- 
tion modes ri^{t) under the translational invariance can be expressed as 

+4A E {^^{tf + / ^k(^)') ^a{t) - HSaO = , 

k^-m'^ + l2X(p^{tf + 12A / rj^r{tf 



(12) 



The derivation and the physical viewpoint in this approach are given in detail in Ref.l). 
Hereafter, let us assume that the chiral condensate points in the sigma direction, that is, 
7^ and = for i = 1 ~ 3. Also, fluctuation modes of the pi directions are identical 
one another, which are denoted as ri^{t), that is, ?7k = = = Vk- When the explicit 
chiral symmetry breaking term, H, is small, the static solutions in Eq.(12) are given as 



lp, = ^o = 4Mf/^\-H/2Ml, r,l' = 1/2^^' + , ryf = 1/2^^' + Ml , (13) 
where the sigma meson mass M^. and the pion mass are defined as 

Ml = -m' + 12A^^ + 12A £ rg^ + 12A j^^ rj^,^ , 



Ml = -m' + 4X^1 + 4A £ 77°,' + 20A £ r]^,^ . 



(14) 



Let us examine the time-dependent solutions around the static configurations. In Eq.(12), 
the time-dependent variables can be expanded around the static solutions of Eq.(13) : 



(15) 



Here, we consider the late time of the chiral phase transition. Thus, Sip{t) is small compared 
with (po. Further, we assume that 1^77^(^)1 -C \r]^\ : 



(16) 



From Eq.(9), it is seen that the two-point function Ga represents quantum fluctuations 
around the mean field. Namely, Ga is small compared with ip^. Thus, from Eq.(ll), we 
conclude that the following relation is satisfied : 



, then , \ipoSip{t)\ > 



(17) 
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Under the above approximations and the small explicit chiral symmetry breaking, the equa- 
tions of motion for Sip{t) and Sri^{t) are obtained from (12) and (15) as 

5^(t) + Ml5^p{t) = , (18) 

5fil{t) + [A{e + Ml) + 24A^o5<^(t)] K(^) = -^AXril^oS^it) , 

5ril{t) + [4(fc2 + Ml) + 8A(^o^99(t)] 5r]J^(t) = -8Ar/^:<^o^<^(t) , (19) 

where 77^ and 5r]^{t) have been rewritten as 77^ ^-^id 577k (^), respectively. The equation (18) 
is easily solved by using the sigma meson mass : 

5^{t) = 5a cos{Mj + a) , (20) 

where 5a and a are determined by the initial conditions. Here, we adopt a = tt for simplicity. 
We further define the dimensionless variables : 

5f]^{T) = ^2714 5rjt{t) , T = Ma/2 , (21) 
/i, = 6X^1/ {k^ + Ml) ■ 5a , K = 2A^2/(k2 ^ . ^^/^^ ^ 

By introducing the above dimensionless variables, the solution (20) and the equations of 
motion (19) are recast into the simple forms as 

7 = 2, (23) 
+ ujI[1- K cos 7t] \5r}^{T) = F^ cos 7T , (24) 




+ujI[1- K cos 7r] j 5r]l{T) = cos 7T . (25) 

If Fo- {Ft^) is negligible, then the equation (24) ((25)) is reduced to Eq.(l). In this case, 
the existence of the unstable solutions for 5fj^{T) may be expected. On the other hand, if 
her (/i-Tr) is negligible or even if the parameters do not offer the unstable regions, the forced 
oscillation may be realized by the effect of F„ {Ft,). 

In the previous paper, the behavior of the time-dependent solutions in Eq.(12) was 
investigated numerically. It was shown that the amplification of the amplitude of the pion 
modes with low momenta was realized even in the late time. However, it was seen that 
the amplitude of the second excited mode was also amplified (See, Fig. 3 in Ref.l)). We 
can now understand their behaviors by using Eqs.(24) and (25) approximately. The model 
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Table I. The values of dimensionless variables based on the parameters used in the previous paper. 
^) Here, S = Sa/ipo 



mode 


UJo- 




ha 


K 


Fa 




n 


= 


4.00 


0.92 


0.766x5 


4.83x5 


12.3x5 


4.09x5 


n 


= 1 


4.08 


1.24 


0.735x5 


2.67x5 


12.1x5 


3.52x5 


n 


= 2 


4.33 


1.89 


0.655x5 


1.14x5 


11.8x5 


2.79x5 


n 


= 3 


4.70 


2.64 


0.554x5 


0.586x5 


11.3x5 


2.41x5 


n 


= 4 


5.19 


3.43 


0.456x5 


0.348x5 


10.8x5 


2.12x5 



parameters m, A and H are adopted so as to reproduce the pion mass, the sigma meson mass 
and the pion decay constant in the static case in Eqs.(13) and (14). The box normahzation 
with the spatial length L {— 10 fm) is applied, in which the periodic boundary conditions 
for the fluctuation modes are imposed. Thus, the allowed values of momenta are restricted 
as follows : kx — {2t: /L)nx and so on, where Ux is integer. The three momentum cutoff is 
also applied to regularize the momentum integration. The fluctuation modes are taken into 
account up to = + + < 8^ corresponding to the three momentum cutoff about 
1 GeV (990 MeV). The values of variables given in Eq.(22) are collected up to n = 4 in 
Table I., where 5 = Sg/^Pq ^ O(10~^) in the late time of the time-evolution of the chiral 
condensate. 

First, let us consider the pion modes. It should be noted that ^ — 2. For the lowest pion 
mode with n — 0, cu^^ — 0.92. Therefore, it would not expected that the amplification due to 
the forced oscillation occurs. We thus investigate the possibility of the parametric resonance 
by setting ~ safely. By noting 2a;,r ~ 2, let us examine the parametric resonance with 
m = 1 in the case Eq.(2). The unstable region is given as 

< e < KuJ^/2 with 7 = 2 = 2c;^ + e , (e fti 0.16) . (26) 

The condition (26) reads together with Table I. as 

5 = Sa/ipo > 0.072 . (27) 

Thus, in the above region for 6a, the parametric resonance occurs and the amplitude of 
the lowest pion mode is amplified. This behavior appears in Fig. 3 in Ref.l). Similarly, let 
us consider the pion mode with n = 1, the first excited mode. As is similar to the case 
n = 0, it is not necessary to take care of the amplification due to the forced oscillation 
because — 1.24. For m—1, the same condition as (26) should be satisfied to realize the 
parametric resonance, except for the value of e : e = 2(1 — uJt^) Ki —0.48. Thus, the unstable 
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region is rewritten as 



6 = 6a/^o > 0.29 . 



(28) 



The parametric resonance is then reaUzed. For m — 2, the unstable region is given as 

-(5/24) • < e < (1/24) • with 7 = 2 = a;^ + e {e^O.76). (29) 

However, (1/24) • hi^LOT^ 0.37(5(7/(/?o)^ is obtained in this parameterization. Therefore, the 
parametric resonance with m = 2 does not occur. 

Secondly, let us investigate the pion mode with n — 2, the second excited mode. If the 
parametric amplification occur, the case m = 2 is realized because 00^^ — 1.89. This condition 
is given in Eq.(29) except for e : e 0.11. However, the upper hmit (1/24) • /i^o;^ gives a 
severe hmit : 0.10 • {Sa/(poy. Thus, the condition (29) can not be satisfied. On the other 
hand, it is possible that the amplitude of the pion mode with n = 2 is amplified by the 
mechanism of the forced oscillation. The value of co^^ {— 1.89) is near 7 (=2). If we neglect 
/i^, the equation (25) is same as (4) and the solution is given in Eq.(6). Thus, we obtain 



Thus, the amplitude increases linearly in time. This behavior is seen in Fig. 3 in Ref.l). 

Finally, let us consider the sigma mode. The amplification due to the forced oscillation is 
impossible because the frequency uj^ is always greater than 2 (=7). Here, let us investigate 
the possibility of the parametric amplification. The lowest mode with n = is the most 
powerful candidate. Since ujcr = 4, the case m = 4 is applicable : 7 = 2 = uJcr/2 + e, where 
e = 0. The unstable region is given as —0{h^) < e < Oih!^) for m = 4 case. Since e = 0, 
the parametric amplification occurs. However, the unstable solution is proportional to e*^, 
in which s is of the order of /?.^ for m = 4 case. From Table L, /4 is about 3.4 x 10"^ with 
5a/ipQ = 10~^. This is very small value. Thus, the amplification occurs slowly in time. In 
the time region under consideration, the amplification may not be seen explicitly. In the 
sigma modes with higher momenta, the parametric resonance would not occur because the 
unstable window is very narrow. 

It should be noted here that the mechanism amplifying the amplitudes of fiuctuation 
modes, i.e., the parametric resonance or the forced oscillation, depends on the model pa- 
rameters. However, in conclusion, it should be pointed out that there is a possibility of the 
amplification due to the forced oscillation as well as the parametric resonance. In this paper, 
the parametric resonance and the forced oscillation are treated separately in Eqs.(24) and 
(25). It is interesting to investigate both the mechanisms in the unified way. The important 
ingredients may be the properties of Mathieu's equation and the Mathieu function. 



(30) 
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